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Abstract — We introduce three types of partial fractional 
operators of variable order. An integration by parts formula for 
partial fractional integrals of variable order and an extension 
of Green's theorem are proved. These results allow us to obtain 
a fractional Euler-Lagrange necessary optimality condition for 
variable order two-dimensional fractional variational problems. 

Index Terms — Variable order fractional calculus, fractional 
calculus of variations, Green's theorem, optimality conditions. 

I. INTRODUCTION 

Fractional variational calculus is a mathematical discipline 
that consists in extremizing (minimizing or maximizing) 
functionals whose Lagrangians contain fractional integrals 
and derivatives. For the first link between calculus of vari- 
ations and fractional calculus we should look back to the 
XlXth century. In 1823, Niels Heinrik Abel considered the 
problem (Abel's mechanical problem) of finding a curve, 
lying in a vertical plane, for which the time taken by a 
material point sliding without friction from the highest point 
to the lowest one, is destined function of height [1]. Abel's 
mechanical problem is a generalization of the tautochrone 
problem, which is part of the calculus of variations (and 
optimal control). Despite of this early example, fractional 
variational calculus became a research field only in the XXth 
century. The subject was initiated in 1996-1997 by Riewe, 
who derived Euler-Lagrange fractional differential equations 
and showed how non-conservative systems in mechanics 
can be described using fractional derivatives [32], [33]. 
Nowadays, the fractional calculus of variations and fractional 
optimal control are strongly developed (see, e.g., [2], [4], [5], 
[11], [14], [15], [16], [19], [22], [23], [24], [25]). For the state 
of the art, we refer the reader to the recent book [21]. 

In 1993, Samko and Ross investigated integrals and deriva- 
tives not of a constant but of variable order [34], [35], 
[37]. Afterwards, several pure mathematical and applica- 
tional papers contributed to the theory of variable order 
fractional calculus (see, e.g., [6], [10], [13], [20], [30], [31]). 
Here, our primary goal is to study problems of the calculus 
of variations with functionals given by two-dimensional 
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definite integrals involving partial derivatives of variable 
fractional order. It should be mentioned that most results 
in fractional variational calculus are for single time, and 
that the literature regarding the multidimensional case is 
scarce: in [3] a fractional theory of the calculus of variations 
for multiple integrals is developed for Riemann-Liouville 
fractional derivatives and integrals in the sense of Jumarie; 
in [12] a Lagrangian structure for the Stokes equation, the 
fractional wave equation, the diffusion or fractional diffusion 
equations, are obtained using a fractional embedding theory; 
and in [28] fractional isoperimetric problems of calculus 
of variations with double integrals are considered. Here we 
develop a more general fractional theory of the calculus of 
variations for multiple integrals, where the fractional order 
is not a constant but a function. 

The article is organized as follows. In Section HI] we 
give the definitions and basic properties of both ordinary 
and partial integrals and derivatives of variable fractional 
order. An extension of Green's theorem, to the variable 
fractional order, is then obtained in Section [Hi] Section [TV] 
gives the proof of a necessary optimality condition for the 
two-dimensional fundamental problem of the calculus of 
variations. We finish with Section [V] of conclusions. 

II. VARIABLE ORDER FRACTIONAL OPERATORS 

In this section we introduce the notions of ordinary and 
partial fractional operators of variable order. Along the text 
L\ denotes the class of Lebesgue integrable functions, AC 
the class of absolutely continuous functions, and by djF we 
understand the partial derivative of a certain function F with 
respect to its zth argument. 

Definition 1: Let < a(t,%) < 1 for all f,T G [a,b], f 6 
L\[a,b], and T be the Gamma function, i.e., 



r( z ) = f e-'i 
Jo 



t z ~ x dt. 



Then, 



t 

«^ (v) /(0= / wj f ^ (f-?)" (f " M /(T)^ (t>a) 
J r(a(f,T)j 

a 

is called the left Riemann-Liouville integral of variable 
fractional order a (-,-), while 



r£ (v) /(0 



1 



r(a(T,0) 



(t-0 



a(x,t)— 1 



f(x)dx (t < . 



denotes the right Riemann-Liouville integral of variable 
fractional order a (■,■). 
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Example 2 ([37]): Let a(t,T) — a(t) be a function de- 
pending only on variable t, - < a(t) < 1 for all t £ [a,b] and 
a certain ngN greater or equal than two, and y > — 1. Then, 

" ' 1 j " r(y+ o(0 + l) ' 
Definition 3: Let < a(f,r) < 1 for all t,T G If 

a i' AC[a, £>], then the left Riemann-Liouville deriva- 

tive of variable fractional order Ot(v) is defined by 



1 



dfi r(l-a(/,T)) 



(t-T)-"^)/^ ( t >a), 



while 



f f >(■'■) f(t) = [ ! f T _ f A«,(v,)-i 

fi > f(t) -J r(a,(T ;f ,)) (T flj 

x /(/] , . . . , f/ i , T,f !+ i , . . .,t„)dx (ti < h) 

denotes the right Riemann-Liouville partial integral of vari- 
able fractional order Ot { -(-, •) with respect to variable ?,-. 

Definition 6: If a .{" J!( ' , ' ) / 6 C'(A,), then the left 
Riemann-Liouville partial derivative of variable fractional 
order OCj(v), with respect to the z'th variable is given by 

«Bp M /(t) = |- ai 4" i(v) /(t) 



while the right Riemann-Liouville derivative of variable 
fractional order «(-,) is defined for functions / such that 

tll" a{r) f eAC[a,b] by 



.( T _ r )-°(V)/(T)JT (*<Z>)- 



d*y r(i-a(T,f)) 

r 

Definition 4: Let < a(f , t) < 1 for all t,T€ [a,b]. If f E 
AC[a,b], then the left Caputo derivative of variable fractional 
order a(-, •) is defined by 



__i_ (t _ t) -«(-,^ /w , t „>,,,, 



while the right Caputo derivative of variable fractional order 
Oc(-,-) is given by 



I 



Let A„ = [a\,b\] x • • • x [«„,£>„], n e N, be a subset 
t = (t\ , . . . ,f„) e A„, and «;(•,•) : [«;,£>,'] x [«,-,£>,■] — > R be such 
that < a; (f; , t) < 1 for all f,-, t G [a;, &,•], i = 1 , . . . ,n. Partial 
integrals and derivatives of variable fractional order are a 
natural generalization of the corresponding one-dimensional 
variable order fractional integrals and derivatives. 

Definition 5: Let function / = f(t\ ,...,?„) be continuous 
on the set A„. The left Riemann-Liouville partial integral 
of variable fractional order «,■(•,•)> with respect to the /th 
variable ?,-, is given by 



1 



(Xi{tt,t)-\ 



r(«*fo,T)) 

X/(*l,---^i-l ! T,?i + i,...,?„VT (f; >fl,), 



1 



_ _->-0((f f ,T) 



r(i t)) ( '' T) 

x/(fl,---,f,-l,T,f, + i,...,?„)c/T (f; >fl,), 

while the right Riemann-Liouville partial derivative of vari- 
able fractional order a, ;(•,-), with respect to the /th variable 
f,-, is defined for functions / such that tj I l b ' ' / G C l {K n ) 
by 



1 /'(t) = -|:r4r ai(v) /(t) 

A/ zi ( 

5rJ r(i -«,■(*, f,-)r 



T-f, 



X/(fi,...,f I -_i,T,^ + l,...,*„)dT (ti<bi). 

Definition 7: Let / e C (A„). The left Caputo partial 
derivative of variable fractional order 0£; (•,•), with respect 
to the /th variable f,, is defined by 

^w /(t)=ai 4-« i (v)|. /(t) 



/ 



1 



r(i-a/(/i,T)) 



-a,(r,-.T) 



X —/(?!,... ,f ; _i,T,f i+ i,...,f„)rfT (f; >fl/), 

while the right Caputo partial derivative of variable fractional 
order OC;(v), with respect to the /th variable f,, is given by 



^■(•••V(t)=-,t ai( -'' ) |-/(t) 



7 m-aKvO) 1 ,j 



r(i - ok{x,ti)) 
•i 

X faf(h,. ..,U-l, T,ti +1 ,...,t H )dT {ti<bj). 

Remark 8: In Definitions |5j [6] and |7j all the variables ex- 
cept f, are kept fixed. That choice of fixed values determines 
a function ft l ,...j i _Yf i+u ...f n ■ — of one variable f,-, 

/fi,...,fj_i,...,^ + i t„(ti) — f(h,- ■ ■ ,ti-i,ti,tj + i, . . . ,t„). 
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By Definitions [TJ g] S |5j |6] and|7] we have 



= ai^tj f(tl,... 7 ti-\ 7 tj,ti+\ ,...,?,,), 
f,^' ' V'l (,-_i,f( + i,...A 



a,A,' ' /fi ',-iA+i fa ( f <) 

v*(v) 



a/A/ ' /( f l:---j f i-lj f !: f i'+lj---i f n)i 



' V*i,...A-_i 1 r i+ i>-A 1 ( f 



a, A,' ' ftu-Ji-lfy+1 t„{ti) 



Cjj/Xi . f^^.^ti-ijtijti+i,.. . ,f„), 



?A>^ ' Vfl fj-lA+l (ft) 

_ C£)^'( ' '/(f[ , . . . ,fi_l,fj,?/_|_l , . . . ,?„). 

Thus, similarly to the integer order case, computation of 
partial derivatives of variable fractional order is reduced 
to the computation of one-variable derivatives of variable 
fractional order. 

Remark 9: If «,(•,•) is a constant function, then the par- 
tial operators of variable fractional order are reduced to 
corresponding partial integrals and derivatives of constant 
order. For more information on the classical fractional partial 
operators of constant order, we refer to [18], [29], [36]. 

III. GREEN'S THEOREM FOR VARIABLE ORDER 
FRACTIONAL OPERATORS 

Green's theorem is useful in many fields of mathematics, 
physics, engineering, and fractional calculus [27]. We begin 
by proving a two-dimensional integration by parts formula 
for partial integrals of variable fractional order. 

Theorem 10: Let j_ < a,(f;,T) < 1 for all f,-,T € [a/,£j], 
where G N, i = 1,2, are greater or equal than two. If 
y , T7i , T72 G C(A2), then the partial integrals of variable 
fractional order satisfy the following identity: 



&, b 2 



g(t) fll C (v) rji(t)+/(t) a2 i? (v) rj 2 (t) 



dt 2 dt\ 



b\ b 2 



m(t)r,C (, 'Vt) + *72(t), 2 i^''V(t) 



dt 2 dt\ . 



a\ a 2 

Proof: Define 



*i(t,T) 



if T < 1 1 
if X > 1 1 



for all (t,r) e [fii,£i] X [a2,b 2 ] X and 




if T < t 2 
if T > t 2 



for all (t,r) e [a\,b\] x [a 2 ,b 2 ] X [021*2] • Since /,g and 77,-, 
z = 1,2, are continuous functions on A2, they are bounded on 
A2, i.e., there exist positive real numbers C\, C2, C3, C4 > 
such that 

|/(t)|<d, |g(t)|<C 2 , |7]i(t)|<C 3 , |TJ 2 (t)|<C 4 

for all t G A2. Therefore, 



/ Fi(t,T)JT+ / F 2 (t 



61 / ^62 / rh 



a\ \Ja 2 \<i a\ 

h 



( fl _ T )«l(fl,T) 



-1 



r(ai(fi,T)) 

( f2 _ T )«2te,T)-i 



-g(t)T}i(T,? 2 ) 



<t/t 



r(a 2 (f 2 ,T)) 



-/(t)TJ 2 (?l,T) 



dx dt2 dt\ 



1 



(h-r) 



ai(fi,r)-l 



-C1C4 



"2 



r(ai(f l5 T)) 

( f2 _ T )«2to,T)-i 



r(a 2 (f 2 ,T)) 



dx jdt2jdt\. 



Because j < G!;(ft,T) <!,;'= 1,2, 



1) ln(fj - t) > and (r ; - t) *^'^ -1 < 1 for n - T > 1; 

2) while ln(fj - t) < and (t t - x) a ' 1 ''^- 1 < (/; - x) r t~ 
for ti — T < 1 . 



Therefore, 



+C1C4 



"2 



1 



T{a 2 {t 2l x)) 



r(ai(fi,r)) 

1 



^_ T )«l(»l.I) 



-1 



(irj c/?2 J df 1 



< i ai li, 2 r 2C3 U. r(a,( fl ,T)) 



dx 



1 



1 



r(a 1 (f 1 ,T)) 

'2-1 1 



(fl-T)'l dx 



-dx 



+ ClC4 \Ja 2 r(a 2 (f 2 ,T)) 

{t 2 -x)^ l dx \ ^dt 2 ^j 



J, 2 -i T(a 2 (t 2l x)) 
Moreover, by inequality 

r(x+i) > 



dt\. 



x z + l 
x+l ' 



valid for x g [0, 1] (see [17]), and the property 

r(x+l)=^r(jc) 
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of the Gamma function, one has 



T{a x {t u x)) 



dx 



' (fi-T)£ _1 </T 



fI _i r(ai(?i,r)) 



+C1C4 



W fl2 r(a 2 (f 2 ,T)) 



''2 



1 



, 2 -l r(0C2(? 2 ,T)) 
&1 / /-i 2 



(f 2 -T)'2 c/T t/f 2 \ dt\ 



1 1 / 1 c 2 c 3 1 r' la ' ( ^ T)+ai(fi ' T) ^ 



+ClC J ^M+^M dz 

\Ja 2 0t£(t 2 ,X) + l 



'2 a 2 2 (f 2 ,T) + a 2 (f 2 ,T) 

f 2 -l 0^(f2,T) + l 



(f 2 -T)'2 a?T u/f 2 \ dt\ 



< / / C 2 C 3 / </T+ / (fi-T)'i AT 
Jaj Wa 2 \ \ 1 Jt\ — \ j 



-1 



/ rt 2 -\ rt 2 _ 

+ C\CAJ dT + J t ( t 2~ x ) h l dT\\dt 2 \dti 



< {b 2 -a 2 )(b\ -a{) 



I b\-\-a\ 
C 2 C3 — 1— ai + h 



+ C1C4 ^1^1 _ 1 _ a 2 + l 2 ^j 



Hence, we can use Fubini's theorem to change the order of 
integration: 



b[ i-b 2 



g{h,t 2 ) ai ff W) r]i{h,t 2 ) 



f{h,t 2 ) al % { ->- ] J] 2 {h,t 2 ) 



1 a\ J a 2 



g(tuh) 



dt 2 dt\ 
1 



a , V{ai{h,x)) 

(fx-T)^)- 1 ^,^ 

rh 1 
+ f{h,h) / 

J an 



r(a 2 (f 2 ,r)) 



J?l(T,f 2 ) 



dt 2 dt\ 
1 



r(«i(fi,T)) 



(*1-T) 



ai(fi,t)-l 



g(h,h)dti 



dt 2 dx 



nb 2 
--2 



172(^1, t) 



1 



t r(a 2 (f 2 ,r)) 



{t 2 -x) a ^ x) - l f{t u t 2 )dt 2 dxdh 



T7i(T,f 2 ) T ^ l( ''' ) g(T,f 2 )rff 2 t/T 



Ja\ Ja 2 



We are now in conditions to state and prove the Green 
theorem for derivatives of variable fractional order. 

Theorem 11: Let < a,-(f,-, t) < 1 - j: for all f,-, T € [«,,£,-], 
where G N, i = 1,2, are greater or equal than two. If 
/,£,TJ G C 1 (A 2 ) and tjl^'^^ll;^/ G C 1 (A 2 ), then 



the following formula holds 

p&l fb 2 p 



'«! Ja 2 
,-b 



^(t)SA? (v) r/(t)+/(t)^ 2 A? ( 



'rj(t) 



r/(t) 



#1 >/^2 



r?(t) 



ir* (v) «(t)*»- 



dt 2 dt\ 
dt 2 dt\ 



'1 % 



rl-<*2(v) 



t"* l ''7(t)<ifi 



Proof: By definition of Caputo partial derivative of 
variable fractional order, Theorem [TO] and the standard 
Green's theorem, one has 



rb\ rb 2 
I u\ J a 2 



g(t)^A? (v) r/(t)+/(t)^A? (v) r/(t) 



dt 2 dt\ 



nbi 
2 



+ ma2 ll^± i{ t) 



dt 2 dt\ 



h fin 



a\ J a 2 



— n(t) T l -^>-) 



vJ /(t) 



dt 2 dt\ 



b\ fb 2 



dt y bl m 



dA 2 
h, 



r?(t) 



j 1 



l-«i(v) 



dt 2 dt\ 



g(t)dt 2 - t2 t b - a ^ ] f{i)d h 



tl D^g(t)+ t2 D^f(t) 



n*b 2 

''2 U b 2 



dt 2 dt\ 



3A 2 



r?(t) h ii; a>W) g (t)dt 2 - t2 ii;^->f(t)d h 



'1 1 b l 



l-a 2 ( v ). 



IV. VARIABLE ORDER FRACTIONAL CALCULUS 
OF VARIATIONS FOR DOUBLE INTEGRALS 
Let (Xi{ti,x), i— 1,2, satisfy the assumptions of Theo- 
rem [TT] We consider the following problem: 
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Problem 12: Find a function u = u(t) for which the frac- 
tional variational functional 



b 2 



subject to the boundary condition 



(tXt)^ A? ( '''Vt)4A? ( '''W)) *2<fti 



(i) 



where y : dA 2 —> M is a given function, attains an extremum. 



We assume that LeC 1 (A2 : 



|; (4 <5, + 2L is continu- 



ously differentiable, has continuously differentiable integral 
t-ll ' , and continuous derivative tj D"'} ' \ i= 1,2. For 
simplicity of notation, we introduce the following operator: 

{«,«!, « 2 }(t) := (t ) «(t) )fl c 1 ^ ( -'' ) «(t)4D^-'' ) «(t) 

A typical example for the cost functional ^ appears when 
one considers the shape of a string during the course of the 
vibration (cf. [3]): 



dxdt, 



where T is the constant tension and cr(x) is the string density. 

Definition 13: A continuously differentiable function u is 
said to be admissible for Problem [T2l if ^D?' ' exist and 
are continuous on the rectangle A2, i = 1,2, and m satisfies 
the boundary condition (|T). 

Theorem 14: If k is a solution to Problem Q~2] then m 
satisfies the Euler-Lagrange equation 



d 2 L {u, oci , «2 } (t) + r, D^ 1 [v) ^ {«, a u a 2 } (t) 

+ ^ (v) 5 4 L{M,a 1 ,a 2 }(t)=0, teA 2 . (2) 
Proof: Suppose that u is an extremizer for J 1 . Consider 
tj e C 1 (A 2 ;R) such that ^.D"'''''^ g C(A 2 ;R), / = 1,2, and 
77 (t)|^ A = 0. We imbed u in the one-parameter family of 
functions {u = u + erj : \e\ < £o,£o > 0}. Define 



J(e) = jf[u 



L(t,fi(t),^ lv, fi(t)4Dg !lv, fi(t))*2*i. 



Then, a necessary condition for m to be an extremizer for 
is given by 



dJ_ 

de 



= 0^ f l [ b2 (d 2 L{u,aua 2 }(t)-'n(t) 

£=Q Ja\ J "2 \ 

+ d 3 L{ u ,a u a 2 }(t)-f n D^' ) r 1 (t) 
+ d A L{u, a u 052} (t) ■ c ai D^ { ''\ (t)\ dt 2 dh = 0. 



By Theorem [TT] and since Tj(t)|,5 Ai =0, one has 



a 3 L{«,ai,«2}(t)^ 1 D^ ( ' ) r7(t) 



- d A L{u, a u a 2 } (t) • ^ 2 D" 2( ' ) T7(t) dt 2 dh 



«1 J«2 



T ? (t)^ 1 D, a i l( ' ) a 3 i{«,ai,a2}(t) 
+ , 2 D^ (v) ^ 4 L{m, ax , a 2 } (t)^j A 2 t/fi . 



Therefore, 



(t) ( d 2 L{u,a u a 2 }(t) 



, 1 D" l( -' , «9 3 i{M,ai,a2}(t) 



-, 2 Dl; { '''d A L{u 1 ai,a 2 }{t) \dt 2 dt x =0. 



Condition (O follows from the fundamental lemma of the 
calculus of variations. ■ 

V. CONCLUSIONS 

Recently, the variable order fractional calculus has pro- 
vided new insights into rich applications in diverse fields 
such as physics, cyber-physical systems, signal processing, 
and mean field games [7], [8], [9], [38]. In this article a 
multidimensional integration by parts formula for partial 
integrals of variable fractional order (Theorem ITOb and a 
Green type theorem with derivatives and integrals of variable 
fractional order (TheoremfTTt are proved. These theorems are 
then used to obtain Euler-Lagrange type equations for the 
minimization of a functional involving derivatives of variable 
fractional order (Theorem IT4l >. Our results generalize the 
recent one-dimensional theory of the fractional calculus of 
variations of variable order [26] to the two-dimensional case, 
i.e., for fractional variational problems with double integrals. 
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